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ABSTRACT 


An experimental and computational investigation has been conducted to determine the off-design 
uninstalled drag characteristics a of a two-dimensional convergent-divergent nozzle designed for a 
supersonic cruise civil transport. The main objective of this investigation was to determine the effects of 
varying nozzle external flap curvature and sidewall boattail angle and curvature on nozzle drag. The 
experimental investigation was conducted in the Langley 16- Foot Transonic Tunnel at Mach numbers from 
0.80 to 1.20 at nozzle pressure ratios up to 9. Three-dimensional simulations of nozzle performance were 
obtained with the computational fluid dynamics code PAB using turbulence closure and nonlinear 
Reynolds stress modeling. The results of this investigation indicate that excellent correlation between 
experimental and predicted results was obtained for the nozzle with a moderate amount of boattail 
curvature. The nozzle with an external flap having a sharp shoulder (no curvature) had the lowest nozzle 
pressure drag. At a Mach number of 1.2, sidewall pressure drag doubled as sidewall boattail angle was 
increased from 4° to 8°. Reducing the height of the sidewall caused large decreases in both the sidewall 
and flap pressure drags. 


INTRODUCTION 

The National Aeronautics and Space 
Administration sponsored a joint High-Speed 
Research Program with US airframe and 
propulsion companies to provide the critical 
high-risk technologies for a Mach 2.4, 300 
passenger aircraft (ref. 1). One task the NASA- 
industry team faced was what type of exhaust 
nozzle would be best for the High Speed Civil 
Transport propulsion system. The technical 
challenges addressed as part of the task 
involved developing a low-noise nozzle with 
high aerodynamic efficiency, low weight and 
long life. Trade studies indicated that a two- 
dimensional convergent-divergent nozzle would 
be best suited to meet the very stringent noise 
requirements of the HSR program. 

The engine cycle chosen for the HSCT 
propulsion system was a mixed-flow turbofan. 
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Although this cycle did have a moderate bypass 
ratio, a mixer-ejector type nozzle would be 
necessary to reduce engine noise. One of the 
nozzle concepts considered is shown in figure 
1 . An ejector would entrain outside free stream 
air at take-off conditions through bypass doors 
integrated with the nozzle. During transonic 
and supersonic cruise conditions, these doors 
would be closed since external air entrainment 
would not be necessary. 

A major problem in designing an 
exhaust system for a supersonic cruise vehicle 
is that its geometric shape must change as 
flight conditions vary. The nozzle exit area ratio 
has to be varied continuously to maintain high 
performance at all flight conditions from takeoff 
to supersonic cruise. Nozzle performance 
cannot be comprised at supersonic conditions 
since payload weight is highly sensitive to 
nozzle efficiency. For example, a 1 -percent 
decrease in nozzle performance has been 
estimated to be equivalent to about an 8- 
percent loss in payload for this class of aircraft 
(ref. 2). 

The nozzle operates with the largest 
area ratio at supersonic cruise. In order to 
reduce flow angularity losses (flow not exiting 
parallel to nozzle center line), the divergent flap 
half-angle and thus, length, are set to provide 
maximum internal performance at the 

supersonic design point for a fully expanded 
area ratio. As seen in the schematic in figure 1 , 
the nozzle would have very long divergent 
flaps. At subsonic cruise, the nozzle would 
operate at a much lower area ratio. 
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Consequently, long external flaps are required 
to minimize the nozzle boattail angle and 
resultant nozzle boattail drag. Another factor 
that can influence boattail drag is the curvature 
at the start of the boattail. 

The nozzle boattail drag can be a 
significant part of the overall drag at subsonic 
speeds. Several studies were also conducted 
that addressed installed nozzle boattail drag 
issues. These studies suggested that nozzle 
boattail drag could be as much as 25 to 40 
percent of the subsonic cruise drag (ref. 3) for 
the HSCT aircraft. Since limited experimental 
data exists for this class of nozzles, an 
experimental and computational program was 
initiated at NASA-Langley to determine the 
uninstalled drag characteristics of this nozzle 
concept. The overall objectives were to: (1) 
determine the effects of nozzle external flap 
curvature and sidewall boattail variations; (2) 
develop an experimental data base for 2D 
nozzles with long divergent flaps and low 
boattail angles and (3) provide data for 
correlating CFD predictions of nozzle boattail 
drag. 

The experimental investigation was 
conducted in the Langley 16-Foot Transonic 
Tunnel at Mach numbers from 0.80 to 1.20 at 
nozzle pressure ratios up to 9. Three- 
dimensional simulations of nozzle performance 
were obtained with the computational fluid 
dynamics code PAB using turbulence closure 
and nonlinear Reynolds stress modeling. 


SYMBOLS AND ABBREVIATIONS 


c D,f 

c D,n 

C D,p 

c D,p,f 

C D,p,s 

C p,f 

c p,s 

CFD 

EXP 

F 

HSCT 


nozzle friction drag coefficient 
nozzle drag coefficient 
nozzle pressure drag coefficient 
flap pressure drag coefficient 
sidewall pressure drag coefficient 
flap pressure coefficient 
sidewall pressure coefficient 

computational fluid dynamics 

experimental 

flap 

height of nozzle, in. 
high-speed civil transport 


2 


Lf length of flap, in. 

L s length of sidewall flap, in. 

M free-stream Mach number 

MS model station, in. 

N nozzle 

NPR nozzle pressure ratio 

rf flap radius of curvature, in 

r s sidewall radius of curvature, in 

S sidewall 

P f flap boattail angle, deg 

P s sidewall boattail angle, deg 


NO ZZ IE D E$IQN§ 

Overall Nozzle Geometry ■ The overall 
nozzle geometry is shown in figure 2. For the 
experimental portion of this study, the nozzle 
was sized to be tested on an existing 
propulsion simulation system in the Langley 16- 
Foot Transonic Tunnel. Thus the nozzle had a 
maximum height of 6.20 inches, a maximum 
width of 6.80 inches and an aspect ratio of 1.10 
(width to height). However, some of the full- 
scale nozzle concepts being considered for the 
HSCT vehicle like that shown in figure 1, had 
aspect ratios of 1.04. Because of this 
difference in nozzle aspect ratio, the height of 
the nozzle was chosen for scaling purposes, 
since this parameter should give the best 
representation of nozzle boattail closure. This 
resulted in the model scale being 8.2-percent. 
From this scale, an appropriate reference area 
can be obtained for subsequent use in 
nondimensionalizing drag in order to produce 
meaningful drag coefficients in terms of airplane 
drag counts. 

All nozzles tested had the same internal 
geometry, because only external nozzle 
geometry parameters were varied. The nozzles 
had a throat height h t of 2.02 inches, an exit 

height h e of 2.70 inches and an area ratio 
Ag/Af of 1.34 (fig. 2). The corresponding nozzle 

design pressure ratio was 5.0. Thus changes in 
performance should only be attributed to 
external flow effects over the nozzle flaps. The 
overall length of the nozzle was 13.14 inches. 

All nozzles were tested with sidewalls. The 
height of the sidewall was fixed to a distance 


American Institute of Aerodynamics and Astronautics 



AIAA 2001-3199 


such that the nozzle internal flaps would not 
unport with the nozzle in the supersonic cruise 
position. In this position the boattail angle of the 
external flaps would be less than 4°. The 
sidewall cross sectional shape and maximum 
thickness were dictated from structural and 
actuation requirements for the full-scale nozzle. 

External Flap and Sidewall Geometry, - 

The parameters used to define the nozzle 
external flap shape are shown in figure 3. 
Nozzle external flap length Lf and flap radius of 
curvature /ywerethe two geometric parameters 

that varied during the test. These two 
parameters were nondimensionalized by the 
nozzle height h m and maximum radius of 

curvature ry max . The baseline nozzle had a 

flap length of 8.42 inches. A shorter flap length 
of 6.74 inches was also tested. Note that the 
overall length of the nozzle, which is 

representative of the full-scale nozzle, was not 
changed when nozzle flap length was varied. 
Nozzles with boattail curvature parameters 
r/rf max from 0 to 1 were tested. A nozzle with 

no curvature would probably be the simplest to 
build since this flap would have a simple hinge 
joint. If curvature were required, some type of a 
sliding mechanism at the hinge joint (fig. 1) 
would be necessary. The term L a j represents 

the portion of the nozzle flap that is a circular 
arc. Nozzle flap FI was considered the 
baseline flap for this investigation. Pretest 
predictions indicated that a nozzle with this flap 
would have the lowest drag. 

The various parameters used to define the 
sidewall external shape are shown in figure 4. 
Sidewall boattail angle (l s and radius of 

curvature r s were the two geometric parameters 

that varied during the test. Sidewall SI was 
considered the baseline sidewall for this 
investigation. 

Nozzle Configurations. Tested . Nozzle 
configurations in which the external flap 
geometry was varied are presented in figure 5. 
The flap radius of curvature parameter was 
varied from 0 to 1.0 for nozzles N1, N2, N3, 
and N4 with the baseline flap length being held 
constant. Boattail angle varied from 11.72° to 
23.44° for these configurations. Flap radius of 
curvature was also varied for the nozzles with 
the shorter flap length as indicated by nozzles 
N13 and N14. All these nozzles were tested 
with sidewall SI. Some tests were conducted 
with a reduced height sidewall with flaps FI 


and F2 to determine whether a drag penalty 
may occur for the full height sidewalls. 

The effect of varying sidewall boattail angle 
and radius of curvature were investigated for 
those nozzle configurations shown in figure 6. 
Sidewall boattail angle was varied from 4° to 8° 
for nozzles N1, N5 and N6 with fs^s.max = 0. 
The sidewall radius of curvature parameter 
r s /r s max was varied from 0 to 0.4 for nozzles 

N5, N8, and N9 with the boattail angle at 6° 
and from 0 to 1.0 for nozzles N6, N6, and N1 0 
with the boattail angle at 8°. 


EXPERIMENTAL PROCEDURES 
Wind Tunnel and Tests 

This investigation was conducted in the 
NASA-Langley 16 Foot Transonic Tunnel which 
is a single-return, continuous-flow, atmospheric 
wind tunnel with a slotted octagonal test 
section and air exchange. The wind tunnel has 
a variable airspeed up to a Mach number of 
1.25. A complete description of this facility and 
operating characteristics can be found in 
reference 4. 

This investigation was conducted at Mach 
numbers from 0.80 to 1.20. Nozzle pressure 
ratio was varied from jet-off up to 9 depending 
on Mach number. Angle of attack was held at 
0°. Reynolds number per foot varied from 

3.8 x 10 6 to 4.2 x 10 6 . 

Air Simulation System. Instrumentation 

The nozzles were tested on an isolated, 
two-dimensional propulsion air-powered 
simulation system that was mounted in the 
wind tunnel by a sting/strut support system as 
shown in figure 7. A photograph showing this 
installation is presented in figure 8. The 
propulsion system is composed of three major 
sections: a nose-forebody, a centerbody 

section, and the nozzle. The nose-forebody 
section up to station 26.50 was nonmetric, that 
is, was not attached to the strain-gauge 
balance. The centerbody section was made up 
of the low-pressure plenum, instrumentation 
section and transition section. The centerbody 
section from station 26.50 to 50.90 was 
essentially rectangular in cross-section and had 
a constant width and height of 6.80 in. and 
6.20 in. respectively. All nozzles connected to 
the centerbody section at model station 50.90. 
The metric portion of the model aft of model 
station 26.50, which was supported by the 
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force balance, consisted of the propulsion 
system, centerbody and nozzle. This simulation 
system is equipped with a flow transfer system 
that is designed to minimize the transfer of axial 
momentum across the force balance. An 
external high-pressure air system provided a 
continuous source of ciean dry air to the model. 
A complete description of this propulsion 
simulation system can be found in reference 5. 

An internal, six-component strain-gauge 
balance measured external aerodynamic forces 
and moments on the model downstream of 
model station 26.50. Mass-flow of the high- 
pressure air supplied to the nozzle was 
determined by a critical-flow venturi system 
located in the air lines external to the wind 
tunnel. Jet total pressure was measured in the 
instrumentation section by means of a 10-probe 
rake. A thermocouple, also located in the 
instrumentation section was used to measure 
jet total temperature. 

The nozzle upper and lower flaps were 
each instrumented with two rows of pressure 
taps with 25 taps per row. On the upper flap, 
one row was located on the flap centerline and 
the other was located outboard near the 
sidewall. On the lower flap, one row was also 
located on the flap centerline and the other was 
located midway between the flap centerline 
and the outer edge of the flap. One of the 
nozzle sidewalls was instrumented with two 
rows of pressure taps with 20 taps per row. 

Data Reduction 

All data for both the model and the wind 
tunnel were recorded simultaneously on the 
facility computer disc. Fifty frames of data, 
taken at the rate of 10 samples/sec for 5 
seconds, were used for each data point. 
Average values of the recorded data were used 
to compute standard force and moment 
coefficients. An appropriately scaled reference 
area for a HSCT transport was used to 
nondimensionalize drag in order to produce a 
meaningful drag coefficient in terms of airplane 
drag counts. 

Faroe Measurements . The balance force 
measurements are corrected for model weight 
tares, balance interactions and flow transfer 
system (bellows) momentum tares. Although 
the bellows was designed to eliminate pressure 
and momentum effects on the balance 
readings, small bellows tares on all balance 
components still exist. These bellows tares 
were determined by testing calibration nozzles 
with known performance over a range of nozzle 
pressure ratios with balance loadings simulating 


the ranges expected for the test nozzles. The 
balance data were then corrected in a manner 
similar to that discussed in references 5 and 6. 

At wind-on conditions, the force balance 
measures all external and internal forces on the 
model downstream of model station 26.50. 
Included are nozzle thrust, nozzle drag (friction 
and pressure drag), centerbody friction drag 
(pressure drag equal to zero since the 
centerbody has no projected area), an axial 
force resulting from a pressure-area term acting 
at the metric break and the bellows momentum 
term. 

Thrust minus nozzle drag was obtained by 
correcting the measured axial force to account 
for the various forces discussed above. Nozzle 
drag was then determined by subtracting thrust 
from the thrust minus nozzle drag term. Thrust 
at forward speeds is determined from the 
measured static thrust for each of the nozzle 
configurations (ref. 6). Nozzle drag includes 
both pressure and viscous forces. 

Pressure Integrations . Nozzle flap and 
sidewall pressure drags were obtained by 
integrating pressures over the respective 
surfaces. Nozzle pressure drag is simply the 
sum of the flap and sidewall pressure drags. 
Note that these terms do not include nozzle 
friction drag which was approximately 0.0001 
for all the nozzle tested. These terms will be 
used subsequently in the analysis of results. 

Data Repeatability 

Data for this investigation was acquired 
within a single wind tunnel entry of the model. 
Short-term repeatability has been quantified in 
terms of a 95% confidence level for nozzle N 1 . 
Examples of short term repeatability of 
integrated pressure drags are shown in figures 
9 and 10 for nozzle N1 flap and sidewall 
pressure drag coefficients respectively at a 
Mach number of 0.90 for seven NPR sweeps. 
Four of the NPR sweeps were performed at the 
beginning of the test and the other 3 were 
taken at the end of the test period. These 
figures show the residuals of the integrated 
pressure drags defined as the difference in the 
individual measured data points from the 
estimated mean of the group of repeated NPR 
sweeps. The estimated mean was the average 
of the grouped data based on piecewise, 3 rd 
order polynomial fits of the individual NPR 
sweeps. Also shown with the dashed lines, are 
the bounds of the 95% confidence interval as a 
function of NPR. In general, the residuals are 
small over the NPR range tested. 
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COMPUTATIONAL FLOW SOLVER 
Governing Equation 

The PAB computer code solves the three- 
dimensional, Reynolds-averaged Navier-Stokes 
(RANS) equations and uses one of several 
turbulence models for closure of the RANS 
equations. The governing equations are written 
in generalized coordinates and in conservative 
form. In an effort to decrease computational 
resources, simplified, thin-layer Navier-Stokes 
equations are implemented into PAB. This 
approximation neglects derivatives in the 
viscous terms streamwise and parallel to the 
surface, since they are typically negligible in 
comparison to the derivatives normal to the 
surface. Extensive details of PAB are found in 
references 7 and 8. 

The flow solver was written with three 
numerical schemes: the flux vector scheme of 
van Leer (ref. 9), the flux difference-splitting 
scheme of Roe (ref. 10), and a modified Roe 
scheme used primarily for space marching 
solutions. These schemes implement the finite 
volume principle to balance the fluxes across 
grid cells and the upwind biased scheme of van 
Leer or Roe to determine fluxes at the cell 
interfaces. Only the inviscid terms of the flux 
vectors are split and the upwind diferenced, 
while the diffusion terms of the Navier-Stokes 
equations are centrally differenced. The details 
and applications of these methods are given in 
references 7 and 8. 

For this study and other typical three- 
dimensional simulations, the solutions are 
computed with the van Leer and Roe schemes. 
An iteration to steady state in a three- 
dimensional computational domain includes a 
forward and backward relaxation sweep in the 
streamwise direction, while implicitly updating 
each cross plane. 

Turbulence Modeling 


model with no damping function was 
implemented in the free stream blocks. 

Boundary Conditions 

The PAB code allows for several boundary 
conditions at the inflow, outflow, free stream, 
and wall and centerline boundaries. Nozzle 
total temperature and total pressure with a 
normal fluid flow angle was used for the jet 
inflow boundary conditions. A jet total 

temperature of 528.67°R was used for all jet 
calculations. Riemann invariants along 
characteristics were used as inflow and free 
stream boundary conditions. A constant 

pressure for subsonic outflow was used far 
downstream as an outflow boundary condition. 
A no-slip adiabatic wall boundary was 
implemented to obtain viscous solutions 

Performance Calculation 

The PAB code contains a performance 
module (ref. 13) that utilizes the momentum 
theorem applied to a user-defined control 
volume to calculate nozzle or aerodynamic 
performance. Quantities such a lift, drag, 
thrust, moments, heat transfer and skin friction 
may be computed for many complex geometric 
configurations and multi-stream flows. Each 
quantity is updated throughout the solution 
development to monitor convergence 

The skin friction force was calculated with 
only the velocity gradients normal to the surface 
contributing to the velocity term of the viscous 
stress tensor. A two point difference was used 
to determine the velocity gradients; one zero- 
magnitude velocity vector at the surface and 
the second at the cell center. Sutherland’s 
formula was used to calculate the viscosity term 
at the surface by extrapolating the static 
temperature at a local cell center to the surface 
and using a reference viscosity and 
temperature condition. 


Turbulence modeling is required to predict Computational Procedure 
solutions for many flow fields. The PAB code 

can perform several turbulence simulations by The computational flow solving code PAB 

implementing either an algebraic or 2-layer, was used to predict solutions at Mach numbers 

linear or nonlinear turbulence models. For this o.90 and 1.20 at a nozzle pressure ratio of 5 
study, the Girimaji 2-equation, nonlinear which was the design pressure ratio. Prediction 

Reynolds stress model (ref. 1 1 ) was chosen 0 f nozzle performance was done for nozzles N 1 

because this model has proven reliable in and N3. 

predicting propulsive flows with mixing, The computational mesh was fully three- 

separated flow regions, and jet shear layers. A dimensional with 3 blocks defining the internal 

modified Jones and Launder form (ref. 12) of nozzle, and 9 blocks representing the free 

the damping function was utilized to treat the stream domain. The far field was located 1 0 

singularity at the wall. A high Reynolds number body lengths upstream and downstream of the 

5 
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aerodynamic nose and nozzle exit, respectively. 
The upper and lower lateral far field was 
located 6 body lengths above and below the 
aerodynamic body. The boundary layer was 
defined for a law-of-the-wall coordinate y+ of 
0.5 on the fine mesh spacing for adequate 
modeling of the boundary layer flow. 

The base grid is a quarter plane 
representation of the experimental model with 
1.57 million grid points in 9 blocks. Using a 
database reduction scheme, a cut 222 grid is 
generated by eliminating every other grid point 
in the i, j, and k directions. This cuts each grid 
dimension by 2, which decreases the grid count 
to 207,437 and substantially reduces the 
memory required to run the flow solver. The cut 
grid can also be sequenced in each direction for 
improved convergence rates and for grid 
assessment. For example, the flow solver uses 
alternating points in the i direction and every 
point in the j and k directions in a 211 
sequence. Generally, a user would begin 
sequencing on the cut 222 grid. A pattern of 
222, 221, 211, and then 111, or no sequencing 
might be used to assess solution behavior as 
more points are utilized in a particular direction. 
The solution is developed until convergence 
requirements are met at each level. Once the 
solution is converged on the cut 222 grid, the 
solution may be extrapolated to the base grid 
and sequencing may again be utilized. 

For example, a converged solution was 
developed within 16.4 hours on a cut grid for 
nozzle N1 at M=0.9 and NPR=5. However, the 
solution was developed to the base level to 
quantify the effect of grid density. Nozzle drag 
decreased a mere 0.00002 in another 30 hours 
of computer time. Since the solution appears 
to be minimally dependent on doubling the grid 
density beyond the cut 222 level, the remaining 
solutions were developed by sequencing on the 
cut grid only. This allowed for quicker solution 
times due to the substantially smaller memory 
requirement. 


DISCUSSION OF RESULTS 
Experimental/Prediction Comparisons 

Nozzle N t A comparison between 
experimental and predicted pressures on flap 
FI for nozzle N1 are presented in figure 11 at 
Mach numbers of 0.90 and 1 .20 at NPR = 5.0. 
A nozzle pressure ratio of 5.0 was the design 
pressure ratio where the flow is expected to be 
fully expanded and peak nozzle internal 


performance usually occurs. The boattail radius 
of curvature parameter was 0.40 for this flap 
and the sidewall was SI. Pressure distributions 
are compared along the flap centerline, and at 
the flap mid station, and at the flap outboard 
station. At M = 0.90, the CFD computations 
tend to over predict the expansion about the 
shoulder of the flap but do accurately predict 
the pressure recovery along the flap. As can be 
seen, the predicted pressures are in excellent 
agreement with experimental data at M = 1 .20 
at NPR = 5. In addition, the pressure 
distributions across the flap are nearly identical 
indicating that nearly uniform flow exists across 
the flap. 

The experimental and predicted nozzle 
pressure drag coefficients are presented in 
figure 12 for nozzle N1 at NPR = 5.0. Shown 
in this figure are the flap, sidewall, and nozzle 
pressure drags where nozzle pressure drag is 
simply the sum of the sidewall and flap 
pressure drags. As would be expected from 
the excellent correlation of pressures, the 
predicted drags are in excellent agreement with 
experimental drag at M = 0.9 and within 0.3 of 
a drag count (one drag count equals 0.0001 
drag coefficient) of experimental drag at 
M= 1.2. 

A comparison of experimental and 

predicted nozzle total drag coefficients is 
presented in figure 13 where the nozzle total 
drag coefficient includes skin friction. Skin 

friction drag coefficient for the prediction is 
computed within the performance module of 
the PAB code. The skin friction drag coefficient 
Cq f is determined using the Frankl-Voishel skin 

friction coefficient as part of the wind tunnel 
standard data reduction system (ref. 6). The 
measured drag coefficient is obtained from the 
force data measurements. As was the case for 
the pressure drag coefficients, excellent 
agreement exits. The total drag coefficients 
agree to within one half drag count (0.00005) at 
M = 0.90 and one drag count (0.0001 ) at 1 .20. 

Nozzle N3 . Experimental and predicted 
pressure distribution comparisons for nozzle N 3 
are shown in figure 14. The boattail radius of 
curvature parameter was 0 for this flap that 
meant it had a sharp corner at the start of the 
boattail. As can be seen, the agreement 
between the experimental and predicted 
pressures is not as good as was the case for 
nozzle N1 at M = 0.90. Consequently, there is 
poor correlation of the flap pressure drag 
(fig. 15) and nozzle total drag coefficient 
(fig. 16). One possible reason for this poor 
agreement may be that the model flap had a 
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small radius at the shoulder. Treatment of the 
corner flow is critical computationally at speeds 
less than Mach number 1. However, at M = 
1.20, where geometrically matching the corner 
is not as critical, agreement between the 
experimental and predicted pressures was 
similar to nozzle N1 . 

Effect of Nozzle Boattail Curvature 

The effects of varying nozzle boattail radius 
of curvature on flap and sidewall pressure drag 
and on nozzle pressure drag are presented in 
figures 17 and 18 respectively. These nozzles 
all had the baseline flap length L/h m = 1.4 and 

were tested with sidewall SI. Results are 
presented at a nozzle pressure ratio of 5.0 that 
was the design pressure ratio. Although 
discussion of results at this pressure ratio would 
generally be applicable to other pressure ratios, 
the relative difference between comparisons 
may vary. Figure 17 also illustrates the typical 
breakdown of pressure drag between the 
nozzle flaps and the sidewalls. Generally, 
varying flap geometry had little or no effect on 
sidewall pressure drag. The breakdown of the 
pressure drags is similar for the other nozzle 
configurations tested. 

The lowest nozzle pressure drag was 
obtained on nozzle N3 with flap F3 that had no 
radius of curvature (fig 18). This was an 
unexpected result because previous experience 
has shown that axisymmetric nozzles with a 
sharp shoulder generally have higher drag 
(ref. 2). However, these nozzles generally had 
shorter external flaps and were designed for 
use with subsonic cruise vehicles. In order to 
try to understand this result, one can examine 
the pressure distributions on these nozzles. 

Pressure distributions along the centerline 
row of the top flap for nozzles N1, N2, N3, and 
N4 are presented at M = 0.90 and 1.20 at 
NPR = 5 in figure 19. Basically what is shown 
is that even though nozzle N3 with no curvature 
had the greatest expansion of flow about the 
nozzle shoulder, it had better pressure recovery 
characteristics than the other three nozzles. 
Similar results were found for the pressure 
distributions along both the mid and outboard 
pressure rows (not shown). 

Figure 20 shows flap pressure drag for the 
two nozzles tested with the smaller flap lengths. 
Except at M = 1 .20, these two nozzles followed 
similar trends as the nozzles with the longer flap 
lengths. At M = 0.90 and 0.95, nozzle pressure 
drag decreased as the radius of curvature 
parameter was increased from 0.1 to 0.4. 


Effect of Nozzle Flap Length 

Varying the length of nozzle the flap is 
illustrated in figure 21 where nozzle pressure 
drag is compared between nozzle N2 and N6. 
Each of these nozzles had a radius of curvature 
parameter of 0.1. As can be seen, the nozzle 
with the longer nozzle flap had lower drag at all 
the Mach numbers tested. This was probably 
due to nozzle N2 having a lower boattail angle. 

Results similar to those shown in the 
previous figure were also obtained for nozzles 
N1 and N5. These nozzles had a curvature 
ratio of 0.4. 

Effect of Sidewall Boattail Angle 

The effects of varying sidewall boattail 
angle on sidewall and flap pressure drag and 
nozzle pressure drag are summarized in figures 
22 and 23. There are no definite trends to 
changing sidewall boattail angle. For example, 
at M = 0.80 and 0.90, nozzle N6 with sidewall 
S3 and p s = 8° had the lowest sidewall 

pressure drag whereas, just the opposite was 
true at M = 1.20. At M = 0.90 (fig. 24), it is 
interesting to note that although sidewall S3 
exhibited greater expansion of the flow about 
the boattail than sidewall SI, it exhibited 
somewhat better recovery such that this 
sidewall had lower pressure drag. Figure 24 
also shows that the flow across the sidewall is 
somewhat nonuniform in nature since 
expansion outboard about the boattail is 
greater for all the sidewalls than at the center of 
the sidewall. However, flow recovers to about 
the same pressure coefficient levels at both 
locations. 

At M = 1.20, sidewall pressure drag for 
nozzle N6 with sidewall S4 was doubled that of 
nozzle N1 with sidewall SI as sidewall drag 
increased from 0.00017 to 0.00035 (fig. 22). 
Examination of the pressure distributions of 
figure 25 at this Mach number reveals markedly 
different flow characteristics across the sidewall. 
Along the center of the sidewall, the flow is 
rather benign and shows relatively little 
differences as p s increased from 4° to 6°. 

Outboard on the sidewall, there is a rapid 
expansion of the flow about the shoulder for 
each of the boattail angles tested. The data 
indicates possible flow separation for the 
sidewalls with 6° and 8° boattail angles. 
Although the full-scale nozzle was design to 
have a sidewall with 4° boattail, it became 
evident as the HSR program progressed, 
sidewalls with boattail angles greater than 4° 
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may be needed for both structural needs and to 
house nozzle actuation hardware. This could 
pose a problem at supersonic cruise where the 
potential now exits for even greater values of 
sidewall pressure drag to occur because of 
sidewall boattail angles greater than 4°. 

Only at M = 0.95, did changing boattail 
angle appreciably affect flap pressure drag 
(fig. 22). Flap pressure drag coefficient 
increased 0.00011 as boattail angle increased 
from 4° to 6°, which was twice as large as the 
increase in sidewall pressure drag. This is a 
somewhat surprising result in that one would 
not expect such a strong effect of varying 
sidewall geometry on flap performance. This is 
because the full height sidewalls tend to act like 
fences and isolate the sidewall flow from the 
flow over the flap. 

Effect of Sidewall Boattail Curvature 

Several sidewalls were tested with varying 
amounts of curvature with both 6° and 8° 
boattail angles. Some typical results for those 
nozzles N6, N8, and N10 with p s = 8° are 

shown in figure 26. Sidewall pressure drag 
decreased by 0.00014 as the r s /r smax was 

increased from 0 to 1 .0 at M = 1 .20. This could 
prove to be significant if sidewall with boattail 
angles greater than 4° are needed as 
discussed previously. Sidewalls with full 
curvature would be feasible for the full-scale 
aircraft since the sidewalls are fixed. 

Similar results (not shown) were obtained 
for those nozzles with p s = 6° except that 

variations in sidewall pressure drag were less 
than for those nozzles with p s = 8°. These 

were nozzles N5, N7, and N9 with sidewalls 
S2, S4 and S4 respectively. 

Effect of Sidewall Height 

The height of the sidewall was fixed so that 
the nozzle internal flaps would not unport in the 
supersonic cruise position. As such, the 
sidewalls are very large unsupported panels 
that account for about 60-percent of the nozzle 
skin friction drag. To assess what drag 
penalties might occur with the full height 
sidewalls, some tests were performed with 
reduced height or cutback sidewalls. The 
height of the cutback sidewall S8 was 
contoured to match the external shape of flap 
FI. Sidewall S8 was then tested with flaps FI 
and F2 to form nozzles Nil and N12 
respectively (fig. 5). 


Sidewall and flap pressure drag and nozzle 
pressure drag coefficients for nozzle N1 and 
Nil are presented in figures 27 and 28 
respectively. Significant reductions in both 
sidewall and flap pressure drag were obtained 
for the cutoff sidewall over the entire Mach 
number test range. For example, at M = 0.90, 
there was a 0.00012 reduction in nozzle 
pressure drag coefficient Cq p and a 0.00023 

reduction in Cq p at M = 1 .20. This would be 

equivalent to at least a 0.00092 reduction in 
drag coefficient for the HSCT vehicle with four 
engines. This does not include a small 
reduction in skin-friction drag. 

The flow characteristics over the nozzle flap 
are different for the two sidewalls. With the 
reduced height sidewall, the flow along the 
sidewall tends to accelerate around the corner 
onto the flap where it can become highly three- 
dimensional. The full height sidewalls tend to 
act like fences in inhibiting any communication 
between the sidewall and flap flows, which 
results in the flow over the flap being more 
uniform. This can be seen in figure 29 where 
pressure distributions along the flap centerline 
and outboard rows are shown for a Mach 
number of 0.90. 

Similar results (not shown) were also 
obtained for nozzle N12; however, the 
reduction in pressure drag was smaller than 
that noted for nozzle Nil. 


CONCLUSIONS 

The results of this investigation indicate the 

following: 

1. Excellent correlation between experimental 
and CFD results were obtained for the 
nozzle with a moderate amount of boattail 
curvature. 

2. The nozzle with an external flap having a 
sharp shoulder (no curvature) had the 
lowest nozzle pressure drag. 

3. At a Mach number of 1.2, sidewall pressure 
drag doubled as sidewall boattail angle was 
increased from 4° to 8°. 

4. Reducing the height of the sidewall caused 
large decreases in both the sidewall and 
flap pressure drags. 
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Subsonic Cruise (M = 0.9) 



Figure 1 . Sketches of full-scale nozzle concept at two flight conditions. 
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Figure 2. Overal nozzle geometry. All dimensions in inches. 
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Nozzles with L/h m - 1.4, sidewall SI 


Nozzle N3, rf/rf max = 0 Nozzles with L/h m =1.1, sidewall SI Nozzles with Lffi m =1.4, sidewall S8 






Nozzle N4, rf/rf max = 10 

Figure 5. Sketches of nozzles with various flaps tested. 
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Nozzle N10, r/r s max = 1.0 


Figure 6. Sketches of nozzles with various sidewalls tested. 
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Figure 9. Short term repeatability for flap pressure drag coefficient, 
nozzle N1, M= 0.90. 
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Figure 10. Short term repeatability for sidewall pressure drag coefficient, 

nozzle N1, M= 0.90. 
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Figure 1 1 . Comparison of experimental and predicted pressure coefficients for Nozzle N1 , 
NPR = 5, L f /h m = 1.4, p f= 16.38°, r f /r fimax = 0.4. 
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Figure 12. Comparison of experimental and predicted pressure drags for Nozzle N1, 
NPR = 5, L f /h m = 1 .4, p f = 16.38°, T f /r lmax = 0.4. 
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Figure 13. Comparison of experimental and predicted total drag for Nozzle N1, 
NPR = 5, Lf/h m = 1.4, p f = 16.38°, r f /r f max = 0.4. 



Figure 14. Comparison of experimental and predicted pressure coefficients for Nozzle N3, < 

NPR = 5, L f /h m = 1 .4, (V = 1 1 .72°, r f /r f max = 0. 
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Figure 15. Comparison of experimental and predicted pressure drags for Nozzle N3, 
NPR = 5, L f /h m = 1 .4, ft- 1 1 -72°, r f /r ffmax = 0. 
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Figure 16. Comparison of experimental and predicted total nozzle drag for Nozzle N3., 
NPR = 5, Lf/h m = 1 .4, p f = 1 1 .72°, r f /^ max = 0. 
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Figure 21 Effect of flap length on nozzle pressure drag, NPR = 5.0. 
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Figure 22. Effect of sidewall boattail angle on flap and sidewall pressure drag, NPR = 5.0. 
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Figure 23. Effect of sidewall boattail angle on nozzle pressure drag, NPR = 5.0. 


20 







Figure 24. Effect of sidewall boattail angle on pressure distributions, M= 0.90, NPR = 5. 
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Figure 25. Effect of sidewall boattail angle on pressure distributions, M = 1 .2, NPR = 5. 
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